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ABSTRACT 


Gamma Ray Dose calculations are done usiing Monte Carlo 
Methods for a spherical shield gecmetry with a straight 
cylinderical duct. All values obtained are for a disc source 
placed at the center of the shield. All the three major 
gamma ray interaction processess namely Photo electric effect^ 
Pair- Product ion and Compton scattering are included in the 
simulation of Photon transport. Results are obtained for 
dose points on the duct axis and also for points off the 
duct axis. The latter shows that the scheme of particles 
scoring used in the present work is not adequate and should 


be modified. 
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CHAPTER- 1 


INTRODUCTION 

Photon Transport studies. ar.e of great importance to 
Nuclear Scientists in shielding analysis and in other allied 
areas such as medicine. In practice^ a few among these shielding 
problems can be handled by mathematical methods. The, then 
obtained analytical solutions or expressions are used to 
estimate the desired quantitites of interest, such as dose, 
in a simple way. 

Often Nuclear Engineers fiace problems which cannot be 
dealt by analytical methods because of irreqular shapes in soxorce 
and shield geometries. Shield geometry depends on the source 
geometry, the specific, use of radio isotope etc. 

Radiation streaming through ducts is also an important 
aspect in the design of shield. This problem was analysed by 
Tsuro et al (11) analytically. The shield geometry considered 
for their study is an infinite slab shield with finite thickness, 
having a cylinderical duct without any bend in it and the 
soxrrce is planar. Apart from research purposes^ geometries with 
openings find extensive use in radiation therapy/non-destructive 
testing. One such geomet-ry has been cho^^jn.. for the present 
study and the geometry is finite in shape. Description and 
the problem of interest, will be described in chap, 2, 
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Though analytical solutions are feasible to some 
extent for duct streaming problems/ which are, in turn/ 
derived with severe assumptions/ one cannot have confidence 
over such solutions unless, they are verified by some other 
methods and ultimately by experimental results. The same 
solutions become inapplicable when source and shield geometries 
are modified. For solving such radiation transport problems/ 
probably Monte Carlo methods are suitable. 

As is well- known that all physical processes including 
the emission of radiation from sources and their subsequent 
transport through material are probabilistic^ in nature/ i.e./ 
one cannot predict with certainity what will happen for 
every individual particle in the prcx^ess, Monte Carlo methcxls 
are used to simulate radiation transport through medim 
mathematically/ through the use Qf various probability density 
functions. Density functions are the basic tools of Monte 
Carlo and are discussed in chap. 2, 

For a. photon radiation transport problem / two broad 
forms of Monte Carlo method are possible (12), 

(i) Ana 1 cxf ue Mont e Car lo t In this method/ the 
schematizsibion followed closely resanbles with the actual 
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physical process. 

In this form/ individual particle/histories are 
followed based on matheimatical probabilities. Particles are 
created by sampling the governing density function of source 
which describes, the behaviour of real particles. Each particle 
is. followed until it survives in the median of interest/ shield 
for e.g.,and then a 'new' particle is started from the source. 
When large number of particles have been accounted/ the 
accumulated data properly handled/ will give the quantity of 
interest, so there is a strong analogy, in this approach, 
between physical and mathematical pcirticles. 

(ii) Non-Analogue Monte Carlo; Non-Analogue method 
is more sophisticated to solve a. problem without any physical 
intution to real particles- or processes. This method attempts 
to solve the transport equation straightaway. 

Although non-aaaJ.ogue method is a versatile approach 
to solve transport problems, it is difficult to form exact 
transport equation for any arbitrary geometry. So more emphasis 
has been shoivn on analogue Monte Carlo method to investigate 
the chosen problem . 
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The main objective of this study is estimating the 
Gamma ray dose rates at seme points of interest, for the 
chosen source and shield geometries. In order to facilitate 
the computations, following simplifications are made in 
developing a Monte Carlo code, for the problem. 

(i) Gamma* rays interact with matter by more than 
ten distinctive processes. However, only three processes 
are significant to Nuclear Engineering problems and other 
processes can be neglected. These are 1 ) Compton scattering , 
2) Photo-electric effect and 3) Pair— production. 

(ii) All scattering interactions aie assumed to be 
isotropic wherever it occinrs in the system. 

(iii) For a history, death can occur when the particle 
becomes absorbed, leaves the geometric region of interest or 
loses its significance owing to other factors (8). 

In the Monte Carlo program, a system built-in function 
random number generator RAN(X) has bedn used throughout this 
work, vjhich generates, random numbers uniformly between 


0 and 1 



CHAFTER-2 


1’4 OiM X Jh ri Ou o 

This. Chapter, deals with the various steps involved 
in the simulation of photon transport. Though there may be 
more than one sampling technique applicable to a perticular 
step, sampling Compton energy, e.g,, here, however, 

relevant sampling methods, say referred by earlier workers, 
are applied for each step of the simulation. Next section 
briefs the description of the problem of interest and the 
necessary steps involved in the simulation are discussed 
in the subsequent sections. 

2,1 DESCRIPTION OF THE PROBLEM 

Present work concentrates on a problem of calculating 
Gamma ray dose rates at sxrnie points of interest for a source 
and shield geometries as shown in Fig 2,1, Source geometry 
is a disc having uniform distribution over it. Shield geometry 
is a solid sphere with a cylinderical duct filled with some 
other material than shield material. Here, main object is 

to compute dose rates aJLong the axis of cylinder and about the 
duct exit also. The latter can be interpreted by computing 
doses at different points on either oidee along the perpendicular 




Fig. 2.2(a) Chosen Co-ordinate System- Spherical Polar Coor- 
dinate (b) Coordinate of a Source Point on the Disc 
Source. 
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to the axis of cylinder , Dose points of interest are shown 
as dots in. the same Fig. on the cylinder axis and on its 
perpendicular also. 

Spherical polar co-ordinate system is chosen to 
represent these, geometries in three dimension and the reference 
axes X, Y and Z are shown in Eig. 2, 2, a. Disc source lies 
in the x-Z plane and cylinder axis coincides with the reference 
Y axis. Shield radius is measured in m.f.p. units which is the 
value at source energy for the shield material. 

2,2 AN OVERVIEW OF ANALOG MONTE CARLO METHCD 

As mentioned earlier/ Gamma photons from a source 
incident upon a shield and thereafter/ its consequences 
in the shield can, however/ be realised mathematically by 
using Analogue Monte Carlo Techniques, Particles coming from 
a source have to be traced individually i.e./from birth to 
death or termination. Further divide the shield medium into 
desired number, of regions or zones according to the choice 
of interest. The procedxire to be followed is shown in 
Pigs, 2, 3. a. and 2,3,b. 

This procedure may be repeated for a large number of 
histories until the desired statistical acurracy is obtained. 




STOP 


Fig. 2. 3. a FLOW CHART FOR MONTE CARLO 
MAIN PROGRAM 







Fig, 2.3 .b FLOW CHART FOR SUBROUTINE HISPEE 
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All steps giv-^n above are discussed in the next few sections 
in relaventfi to the chosen problem, A random process or 
event can be described mathematically by making use of the 
idea of probability density function (pdf ) and cumulative 
distribution function (cdf) which are defined as follows: 

(i) Probability Density Function : p(x) describes the 
relative frequency of occurance of its random variable 
X (12), X can lie anywhere between -cd and + oo,' A pdf has 
the following propertiea: 

1) 0 < p(x) for all X 

+CD 

and 2) p(x) dx = 1 

(ii) Cumulative Distribution Function. ; C(x) givea the 
probability that a random variable x has. a value less than 
or equal to a fixed value. It is defined by 

X 

C(x) = / p(x) dx = probCx^X) ,,,,,(2,1) 

-OD 

cdf haa the following propertiea: 


1 ) 


lim 

X + oo 


C(X) = 1 
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lim 

X C30 


C(X) = 0 


To sample a random value x from a given pdf# first 
select a value C(x) which is ?nerely a number# essentially 
from numbers distributed uniformly between 0 and 1, then X 
can be obtained directly from equation (ij- ) provided the 
corresponding pdf is a simple one# otherwise some other 
methods are to be followed to saifiple X. 


2*3 SOURCE PARAMETERS 

In order, to initiate the problem# one has to start with 
source particles i.e*# photons and their parameters* Because# 
it is essential to know their spatial point of origin and 
direction of motion* A disc source is madeup of number of 
point sources* Since the choice of so-urce parameters for a 
photon is not the same for different source geometries and for 
different source distribution like isotropic, cosine distribution 
etc*# only uniformly distributed disc source and isotropic 
distribution for point source have been assumed throughout 

this. work. 

i) Sf^Tection of Space Co-ordinatea : Consider the problem of 
choosing spatial co-ordinates for source particles. The 
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co-ordinates (x,0,z.) of any point on the disc can be found 

interma of a selected randtm value R betv;een 0 and radius 

of the disc^ and an angle randomly selected between 

- A and + . R and y are shown in Fig, 2.2.b. The 

probability density function describing the distribution 

of sources distributed on a disc as a function of radius 
2 

is 2 r- R/ R^ , where Rj- is radius of the disc^. and hence ^ 
the cumulative distribution fionction from equation (1) 
is (3) (8) 


R 

C(R) = / 
o 



R /R 


2 

I 


using a random nxinber 'r' between 0 and 1 instead of C(R) 
to select a random value of R directly from the above 
eqnation/ it gives 


2 2 
r = R /Rj 

It follows that 



.,*.( 2 . 2 ) 


From the fig, 2.2,b/ it can be seen that for a given 
value of R, the source positions are equally distributed 
in azimuthal angle y between - A snd +7^ • Then selecting 
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for Q random number 'r' , it gives 

r. = / dy/27^= (■f>+/r)/27N 

or Y = ^(2r-l) (2,3) 

Nov/ valuea of R and ^ are known, so point Tix^z) can be 
located on the disc from the equations x = R Cosy^nd 
z = R sin V . 

ii) Selection of Direction Parameters : After aasigniog . 
spatial point of a source particle^ or more precisely, after 
locating a point source on the disc, next consider the 
problem of assigning direction parameters to source particles. 

A disc source is symmetric w,r. to its axes i.e, it emits 
radiations equally on either sides. Further, a point source emits 
radiations into entire 4 /c steradians of solid angle. Here 
again ^one can consider only a half of the sphere by taking 
the advantage of symmetry about a plane. So now the problem 
has reduced to a half sphere pointing towards +Y direction 
and hence, if a point is chosen on it, that point will give 
the necessary direction to the soxirce particle. Of course 
that point may, however be expressed interriis of direction 
cosines using a polar angle g and an aaioivtbal angle 0 • 
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4) and 0 are shown in Fig. 2, 2, a. Selection of ^ and @ can 
be done in the following manner. 

1 ) To select a polan angle ' 0 ' randomly anywhere 

between and 0/ the pdf. interms of cosine of the angle 
'k* - between -1 and +1 - is j dk (3 ) and from this pdf, a 
random number r selected for the cdf gives directly. 



or k = 2r-l = CPs0 (2.4) 

2 ) To choose an azimuthal angle '0' between 0 and a , 
it is assumed that 0 varies evenly in the interval 0 and k • 
So, for a random number 'r', choose to be equal to^r. 

Finally, direction cosines for the soiarce particle 
in spherical polar co-ordinates are obtained from the 
following equations using the values of 0 and 0 , 


U = Sine Cos0 

Uy = Sine Sir0 ...••(2,5) 

U = Cose , 
z 
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dir.ectd.on cosines, along X/Y and Z 
directions respectively. A subroutine is written to compute 
space and direction co-ordinates of soxxrce particles using 
equations (2.2)-(2.5). 

2,4 PATH LENGTH SELECTIOH 

Thusv s s.et of source particle's parameters has been 
selected. Next, the path length from the source to the point 
of interaction has. to be determined. The path length alongwith 
parameters of initial direction cosines defines the point at. 
which interaction occurs. But, two mediums l^ve to be 
considered for the geometry. Let duct be filled with air and 
the remaining part of the shield is made of lead. Only 
homogeneous medium has been considered wherever lead or air 
is present in the geometry. So before sampling the path length 
considering the duct, photons emanating from the source will 
fall into one of the following three cases listed below: 

(i) A photon originating from the source may escape 
through the duct i,e,, photon without even making any 
interaction with shield comesout of the duct freely passing 
through different mean free paths as. shown in Fig, 2. 4, a. 





Fig. 2.4 Possible Particles Line of Flights . 
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This contributes a part to the xoncollided or unscattered 
dose at the source energy* One m.f.p. in air at 1,25 MeV 
is approximately 13x10 cm whereas the whole spherical shield 
radius is approximately 9,0cm only and therefore, one can 
neglect the attenuation due to air throughout the particles 
line of flights. 

Solving equation of a straight line and that of a 
cylinder will give where '^ 2 .' shown in lig, 2, 4, a. 

The co-ordinates of that arbitrary point Q(x' ,y' ,z' ) can 
be calculated using the following equations namely. 


x' = X- + U^kj^ 

y* = y + Uyk^ .,...(2.6) 

Z* = Z + 

If / (x' + (y' 5^ + (z' is greater than the maximum 

radius of the sphercial shield, th-sn that photon will escape 
through the duct, otherwise it will not. 

ii ) The situation is shown in Fig, 2,4 ,b. Photon 
does traverse some distance •}C 2 '-found in the same manner as 
I ^ through the duct and the attenuation in air medium. 
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obviously-/ can be neglected, since 1^2 Hence particles 

CO— ordinatea have to be modified in order to keep track of 
the particle. Substituting k 2 in place of k^ in equation 
(2,6) will yield the latest position of the particle with 
initial energy Eq, from where the particle will travel into 
shield by a distance equal to the selected value of the 
path length. However/ direction parameters are retained for 
this purpose. In some other situations i.e,, when is 
greater than duct radius, it may happen to travel some distance 
in air from its origin to the surface of inner cavity sphere. 
This is not shown in the figure. 

(iii) In this case, photon does traverse no distance 
in air. Of course, this will arise only when both source and 
shield are not separated by any distance from each other 
as shown in Fig. 2,4.c. and ofcourse, the chance of occurance 
of such case is very less unless the geometrical considerations 
are modified. So photon may be allowed to travel into shield 
directly just by sampling the path length. The selection 
procedure for the path length 'S' is as follows. 

The probability of first collision between S and 
the line of flight is given by (2) 


S + dS along 
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p(S) dS = ds 

where ^ Is the total linear, attenuation coefficient of the 
shield materialvAhieb givea the probability of collision 
per unit length. To sample ordf is given by 

c.(S) dS 

Select a value e(s) using a randwi number, r between 0 and 1 
above equation becomes 

r = / ^dS = 1 - 

0 

It follows that 

S = - i In(l-r) 

But 1 - r is distributed in the same manner as r. Therefore^ 

S=:-.ilnr (2,7) 

The selected path length 'S' together with 
direction cosines gives the point at which the photon interacts 
with the shield. Interaction may be any of the three events 
namely Compton scattering. Pair- Production and Photo electric, 
effect. Sampling techniques applied for these three events 
will be discussed later. 
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2.5 NEW DIRECTION PARAMETERS 

This section discuss the problem of finding new 
direction of a photon after it encomtered a scattering or 
collision event. Consider a Ccmpton scattering event occurs 
at aome point in the shield. Once an elastic collision occured 
in a rnddium^ its consequense results change in photons 
direction with loss in its energy. The situation is shown in 
Fig. 2.5 .a. A method of sampling photon energy after collision 
will be discussed in later sections^ . In this section^ 
only particles new direction will be determined for a given 
scattering or deflection angle and its local azimuthial 
angle as shown in Fig. 2.55,aand 2.5;,b* ©^ and 0^ are 
measured w.r. to particles initial direction. The new direction 
parameters v;ill be very useful to maintain the track of 
a particle from its birth w.r. to a particular referereae axes. 

Having known old values of © and 0 and the values of 
© and 0 / the new direction pirameters U ' and U * 

are determined with the help of the following four equations 
where U ' ^ U ' and U' are the modified direction cosines 

w « 111 A * 'y 2 

w.r. to different , axes which are parallel to the 




scattered electron 


Fig. 2.5.3 Conptcm scattering: an elastic collision between 
a photon and electron. The local orthogonal (Cth»?) 

are constructed such that 5 lies along the original direc” 
tion of the photon’s travel. 



CHcte*,-:) Particle's 'local' angles of scattering, 6^ is the 
deflection angle, f is the azisnithal angle. Usually the 
angle • can be assumed to be randomly distributed in the 
range 0 to 2r. 

Figo 2«5ob Initial ai^d Final Direction of Particle 
in Scattering 
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original reference axes. X, Y and Z respectively, 

Cos©^ = Cos0Co^ -► slnesinS^Coa 

Sin©^ = /I - cos©^- Gos©^ 

Cos0j^ = Cos(i2^ ~ 0 )Cos 0- Sxn(0^ - 0)s±n0 

and Sin0j^ = Sin(j2^ ~ 0)Cos0- Cos.i0^ - 0)S±n0 

©^ and 0 ^ are the modified polar and azimuthal angles 

respectively. The above, foxur equations are taken from Wood 

and Ref, (6), According to Ref. (12), these equations are 

obtained by making use of trignometric relations. Once ©^ 

and 0 ^ are known, U U„‘ and U ' are readily found by 
n -jc y z 

substituting 0^^, 0 ^ in place of © andj^in equation (^.5). 

Since the purpose of the above equations is to obtain the 
new direction parameters, much interest is not shown in the 
details of the derivation of them. The derivation may be 
seen in Ref. (12), In addition to the above method, few other 
methods are also proposed in Refis, (3) and (8), 

2.6 ABOUT THE DUCT 

With the knowledge of photons new direction parameters 
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and the new energy E < 1 "" energy), sample 

a path length 'S'« from equation (2.'7 ) and calculate the next 
interaction point using equations given below 

Xn = X ' + U ' ■ a • 

Yi = y* + u* s* 

and z, = 2 ' + U* S' (2,8) 

The next interaction point may face^ apart from the type 
of interaction, any one of the following three situations- 
given below; 

(i) It may lie even outside the sphere or shield 
or so to say, photon escapes from the shield with energy E, 

Of course, this can happen at any time independent of the mjmber 
of collisions. If a photon . escapes from the shield then 
ultimately its life history will be terminated. 

(ii) It may lie within the shield^ Here it is 
essential to decide the type of event to be aceounted, 
depending upon that only, photon may either be allowed to 
survive or not to survive. 
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(iii) Photon may encounter a situation as shown in 
Fig, 2,4d-. i.e, photon crossing the duct. Here sampled 
path length 'S'* kept linaltered but it traverse an additional 
distance of 'k" ' in air along that direction while crossing 
the dxx:t. Whether photon crosses the duct or not can be found 
easily from the values of ^Cj^y^ and So if a photon 
crosses the duct^ then its actual point of interaction will 
be found just by substituting (S' + k" ) instead of S' 
in equation (2,8), 

Further when a photon is allowed to cross, the duct 
the point of interaction may or may not lie within the shield* 
In the former case^ photon just crosses the duct and in the 
later, case photon escapes through the duct; ultimately its 
termination* 


2,7 SAMPLING TECHNIQUES 

2*7,1 CcmptQT^ .qn;=^t-.t.firina t This effect is the scattering of 
a photon by a free electron which is at rest. In this 
interaction.^ photon loses, some of its energy from its 
original value to sane other value E (E^> E) and is 
deflected by an angle 0^ from: its original direction of 
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motion (see Pig, 2, 3, a). The relationship between E and 
©o ©repressed interms of Compton wavelengths units namely 
% = 0.511/ E , where E is the photon energy in MeV^ is 

A— = 1 - COS©^ 

® O 

where and ^ are the wavelengths before and af ter 
scattering respectively and ( X - gives the increase 
in wavelengths which corresponds to the decrease in energy 

In Monte Carlo calculations^ it is recessary to know 
after a Cempton scattering ^ the value of 7\. and ©^. Sampling 
from Klein-Nishina (K-N) scattering probability density 
function expressed interms of X = ?> / Aq v/ould suffice this 
requirement . 

Though there are many techniques available to sample 
X from the K-N pdf. Wood and Blomquist (1) suggest, that the 
sampling technique valid for any incident photon energy and, 
ofcourse, the most widely used also, is that of Khan (7). The 
details of K-N can be seen in He£.(12) and the derivation of 




X Value 


fig. i.6 flow .cy^oxt.. for Kahn's procedure for saa|>ling 
the K-N distribution to find the wavelength of a gama pho- 
ton gfter a Compton scattering. 

In is photon wavelength before scattering (assiaaed known) « 
is photon wavelength after scattering and X*R»Xu 4 .|/Xu. 
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Khan sampling technique may be seen in Ref,(l), The procediore 
for sampling the K— N pdf to find the n&i wavelength ' X ' 
is shown in Rig. 2.6, 

Using the selected value of X shown above sampling 
procedirre^ CosG^ can be obtained from the equation 

CosG^ = 1,0 - X + A 
0*00 

Now the values, of a and 9^ are known. Next, the azimuthal 

angle 0^ may be selected in the following manner.. As is 

obvious, 0Q can have any valuesb between 0 and 360°, see 

Big, 2,5,b, This angle may be obtained by sampling at random, 

assuming 0^ is distributed evenly in that range, as 

0 = 2Ar, where r is a random number between 0 and 1. 

^o * 

2,7,2 Pa lr~ Pr oduct i on : This is an important reaction for 
photon energies greater than 1.022 MeV, In this reaction'’, 
the incident photon disappear a in the field of charged particles 
and an electron - positron pair appeaifewith total energy equal 
to that of photon. The energy (E - 1.022) MeV goes into 
the kinetic energy of electron and positron pair. They 
are attenuated in the medium. The electron eventually 



ttrigiiiAl Y^photvi) 

{£>1.022 (<eV) 


positron combines with\ 
stray electron and is 
annihilated 


0.511 Me 



pair production interaction 
occurs here: positron-electron 
pair created 


electron eventually 
absorbed in material 


0.511 MeV 

2-anni^ilation quanta created in 
random, but opposite, directions 
Fig. d. Scheme of the pair production process. 



(Hote*.-') ihe 2-annihilation ganaaa photons are aesuned to 
saerge ie ©tropically. That is, with respect to the ortho- 
gonal axes (x', y*, s'), the arimuthal angle ♦p is uniform- 
ly distributed in the range 0 to 2w, and the cosine of ©p 
is uniformly distributed in the range -I to I, 

2^7^13 j)j_r0Cti.on of Anniltilsttion Gamma. s Afteir 
Pair- Er oduct i on 


Fig, 
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absorbed while the positron combines with an electron in the 
attenuating material and in this process, two Gamma quanta 
of 0,511 MeV (rest mass of electron) are released. 

In Monte Carlo calculations^ for accounting a 
Pair-Production event, two assumptions are made. Firstly^ 
annihilation Gammas are released at the position at which 
the Pair- Production occurs. Secondly, the two Gammas are 
released at random in mutually opposite directions. They 
are shown in Fig. 2, 7, a. 

Now, the initial direction of annihilation Gammas 
has to be sawmpled. The sampling is done in the following manner* 
Further, it is assxned that the angular distribution of 
Gammas to be isotropic. That is, the cosine of 0p, the polar 
angle, is uniformly distributed in the range -1 to +1 and 
the azimuthal angle 0 ^ is uniformly distributed in the range 
between 0 and 2 k. 0^ and 0 ^ are measured w.r. to parent 
photon's direction and they are shown in Fig. 2.7.b, 

According to Ref. (12), to reduce C.P.U. time, only 
one annihilation Gamma from the pair can be followed; as 
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compensation / its weight must be mxiltiplied by two* Concept 

of weight will be seen in Chap, 3 . The error resulting from 

the above approximation is likely to be very small because 

of the random orientation with which annihilation quanta are 

released. Finally, the value of © and 0 are selected from 

P P 

the equations. 

Cos© = 2r. “ 1,0 

p X 

and - 27?r- 

P ^ 

where end are random numbers. Using these values the 
new direction parameters w.r. to the reference axes, can be 
computed using the equation given in section 2 .5 . 

2,7,3 Photo Electric Effect ; It is an absorptive event and 
in this process, the incident photon is absorbed and 
depending upon its energy, an inner orbital electron (i,e. 
fron k shell ) is ejected. The vacancy created is filled by 
a transition of one of the outer shell electuons accompanying 
either a fluorescent emission or an Augar electron. This effect 
is. very important for low energy photons in materials with 
high atomic nijmbers (Z), in this Monte Carlo calculations for 
Photoelectric, effect it is assumed that all transitions 
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are accompanied by the emission of fluorescent lines only. 
Thus^ Augar electrons are completely neglected throughout this 
work. 


The energies corresponding to K and L shells are 
called edge energies because the Photo- afciectric cross section 
elsewhere decreasing continuously shows a sharp 
discontinuity (edge) at each of them (2). The behavioiar at the 
K and L edges are shown in Fig. 2. 8. a.. 


The energies and compared with the incident 
photon energy serve to determine the probability of K or L 
shell ejections. So when a photon of ener^gy E > Ej^ is absorbed 
there is a likely chance of K-shell electron ejection with 
kinetic energy (E - E^) . Eurther the vacancy created ia 
filled by an electron from any one of the outer shells (l^ 

M or N) making a transition accompanying a fluorescent 
emission. The possible K and L fluorescent anission or lines. 
Ktx , K« . Kfl and F, are shown in Fig. 2.8 .b and their 

energies are given in Table 2,1. For elements with 
L - fluorescent emission is also possible. 


If the incident energy E is greater than then 
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either K or L shell ejection is possible. Qfcourse outer shell 
ejection is aiaso allowed. Now, one has to decide which ejection 
is to be considered. So some three nmbers and 0^, 

which will determine the relative probabilities of K, L and 
outer shell ejections for E v.E^ and 0^ also define the 
chances of L or outer, shell ejections, for E^< E< E^,, are 
required. 

The values of 0 , 0 and 0^ are taken from Ref, (5) 

o ij Jx 

and are shown in Table 2.1, These Vcilues are computed 

using the Photo electric cross sections at various edges. 

Further, the yield for a shell i is the total probability 

of emission accompanying electron transition from all outer 

shells to a vacancy in shell i. The individual yields for K - 

lines and L - lines are shown in the same table, A random 

tin 

number lying between implies i shell electron 

ejection and similarly a random number between 

and Yj_/ 0^ implies i^^ fluorescent emission. Here also, 

the particle’s statistical weight is modified and this concept 

will be discussed in Ghap. 3. Next, the photon with new 

fluorescent energy has to be followed, Thrs, as done before 

in case of Pair- Production, those two angles are sampled in 
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Table 2.1 Final values of j2. and F. 

I --- r- - - - - — - - -1^ 


i 

Shell 

Edge 
Energies 
in MeV 


Y. 

1 

Flu. Ene 
in MeV 

Flu, 

Lines 

1 

0 

0,014699 

0.242 

0,0 

0.0 

0 

2 

L 

0.014699 

1.00 

0,3032 

0.012217 


3. 

K 

0.088004 

2.854 

2.074 

0.074969 


4. 

K 

0,088004 

3,957 

3.127 

0.072804 

^<2 

5 

K 

0.088004 

4.570 

3.712 

0.084777 

Si 

6 

K 

0.088004 

4.715 

3.850 

0,087320 

K,2 


Z = 82 



the same manner and finally# as usual# the new direction 
parameters can be computed using equations given in section 2,5 


Although# the method to be followed for a Photo 
electric event is described clearly in Sefa, <4i and (5) 
the derivation and the method of obtaining j3^s ' and s * 


are 


not explained. 



CHAPrER-3 


RESULTS 

In this Chapter, resiilta, obtained by making use 
of the different sampling techniques that are to be followed 
at various, stages of photon transport simulation seen in Chap. 2 
have been presented. All necessary data, are of Israel (10) 
and that of mass absorption coefficient of air or tissue 
are from! Ref, (12). in the development of a computer, code 
for the chosen problem, three • cases have been considered 
to observe the improvement of results of one case over 
the other and of course, to attain the closeness of exact 
physical situation also, 

Observatiorshave been made firstly, by forcing all 
interactions to be Compton scattering (exam), secondly, by 
considering only two types of interactions namely Compton 
scattering and Pair** Product ion event (com+ppd), and thirdly 
by including all the three events viz,, Compton scattering, 
Pair-Production and Photo electric effect (com+ppd+pee). 

Next section gives a brief idea about the concept of survival 


weights 
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a.l CONCEPT OF SURVIVAL WEIGHTS ^ 

In Monte Carlo computations, this idea, of statistical 
v;eight has been included to account the absorption effect 
after a scattering, A particle starts from the source with a 
weight of unity. Obviously it has not encountered any 
interaction in the medium so far. The particle weight after 
a collision is the weight before collisiooj multiplied by 
the survival ratio Pc/^ which is the ration of 

macroscopic scattering cross section to the macroscopic total 
attenuation cross section. 

Thus after each scattering event the weight is 
adjusted in the same way to account the effect of absorption 
of photons in the scattering. Here, one can choose R ussian 
Roullette method (8) to decide whether the particle survives 
or not , ^ 

3*2 COMPUTING DOSE 

As mentioned earlier, the five mfps' shield is 
divided into five regions. Dose conputations have been 
made for. each mfp or region at points on the cylinder axis 
and along the circumference of the circle with radius equal 
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to that of spherical shield in the Y-z plane to accoiint for 
the dose pointa along the perpendiciolar- axis of the duct. 

The right half circle is divided into 18 equal parts w.r, to 
the origin in step of 10 degrees each and hence doses computed 
in each interval corresponds to seme amount of height in the 
perpendicular axis. 

Though results are obtained for all the 18 intervals / 
assuning symmetry one can limit to only either upper half or 
lower half of the semicircle. From the observations^ it .is 
foxind that there is not significant difference in values betwesen 
the two halves. 


VJhen a particle crosses a mfp,the associated particle 

2 

current is converted into flux units^r say T (no, /cm - sec) 
and then dose rates in miV'hr can be calculated using the 
relation namely. 


(E) 

Dose Rate = 0.05767 X E X X Tl 

where 0.05767 is the constant which converts the unit of 

u_ 

E --- T into dose rate mR/hr units, E is the En'ergy of the 

f P 

photon in MeV and --- is the mass absorption coefficient 

^ 2 

of air or tissue at energy E in cm /gm. 
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3,3 MONTE CARLO OUTPUT 

The three cases mentioned above are given 
below along with the results obtained in each case, 

(i) Here all interactions are forced to be 
Comptons scatterings i.e,# absorptive events cannot occur 
directly. Thus the particle weight is adjusted for by a 
factor of survival ratio after, each collision. Results 
obtained are shown in Pigs, 3,1 and 3,2, Fig. 3,1 shows the 
distribution of doses rates at each mfp for various angles 
i.e., along the perimeter of the circle. 

(ii) In thiS/ all interactions are forced to be 
either. Compton scattering or Pair— Production event* Cf course 
if the energy of the photon is less than 1,022 MeV then all 
interactions are ultimately forced to be Compton scatterings. 
If E 1,022 MeV^ then I^ir-Production is considered to occur 
if r where r is a randcsm nun her between 

0 and 1 t otherwise it will not occxar and here the survival 
ratio is given by 

Survival ratio = 

As describes in (i) above^ results are shown in 




1 2 3 4 5 

Shield Thickness ( mfps') 


Fig. 3.1 Variation of Dose with Mean Free Path (Com) 



Dose Rate, (mR/hr) 



Fig.3.2 Variation of Dose with Angle ( Com ) 
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Pigs. 3.3 and 3,4, Conparision study will be done later, 

(iii) In this case ©hoto Electric, effect is also 
included in the simulation along with other two process 
namely Compton scattering and Pair- Product ion. Now one has 
to decide the type of interaction to be allowed to occur 
at the interaction point at seme energy, say E. This can be 

done by comparing a random number 'r' with the ratios 
FpEl/yT* '^PPD/^T, ^ lies between Ppgg/ 

and p^, (Say /^ppj^ 4 < JJq at E), then the Compton 
scattering will occur. Likewise the choice of the event may 
be decided at any energy E. 

Here also the particle weight after an event is 
modified by the factor survival ratio + U^)/ pL. In 

the case of Photo electric, eftect, if fluorescent emission is 
allov/ed then the particle will survive with that fluorescent 
energy Otherwise its life will be terminated. Results obtained 
in the case are showrr in Pigs, 3, 5 and 3.6, 

3.4 COMPARISION OF MONPS CARLO OUTPUT WITH ANALYTICAL METHODS 

A semi— analytical approach for radiation streaming 
through a cylinderical duct has been made satisfactorily by 
Tstoto et.al, (11). Eirst this method was applied to the case 



Dose Rate , ( mR/hr) 



Shield Thickness (mfps*) 

Fig. 3.3 Variation of Dose with Mean Free Path (Com ♦PPD) 
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10 ^ 
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Fig. 3.4 Variation of 




Shield Thickness (mfps’) 


Fig. 3.5 Variation of Dose with Mean Free 
Path(Com+PPD*PEE) 


Dose Rate , (mR/hr ) 



Fig.3£ Variation of Dose with Angle (Com^PPD ♦PEE) 
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of Gamma ray streaming after single scattering in a duct for 
a. disc aource at its entrance. Later , it waa modified to 
include multiple scatterings also. The dose due to single 
scatterings or singly scattered dose at the duct exit in 
given by 


SCA 


^ L +00 

■■TkR^ / dz/ --- 


0 


Po 


d £__ __1 

dfi- Dq 


exp(-^do +>^1^1) % d£ 


where. L and Rq are the length and radius of the duct 
respectively 

PqS initial attenuation coeficient of photon 

attenuation coefficient of photon after scattering 
dQ: sampled path length 

dj^: distance between scattering point and detector 
minus the path across the duct 
DqJ distance between the point of origin and the 
scattering point 

distance between the scattering point and the 
detector 

2 

f : an integral parameter and r = Rq + f / 
the piboton energy after scattering 
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is the energy abs-orption coefficiert of air at 

J 

energy 

The total dose at any point on the duct axis, is the 
s^Jm of unscattered dose and scattered dose ^ and the unscattered 
dose at any point 'L' units, away from the disc, source is 
given by (9) 

(Eq) ln(l+R^L^) ^...•(3.2) 

X 0.0567 

where is the flxox and is the initial energy of photon 
( at source ) . 

Estimation of dose rates by analogue method is as 
follows. Attempts were made firsts to calculate the doses 
along the duct axis and later for points off the axis. As 
such, no point detector calculation^ say a particular point in 
the shield, can be done straight away Using this method 
but the normal procedure is to consider some distance +4id 
about that point along the boundary so that the sm of the 
individual contributions of each particle's dose within that 
total 2 2id length will, be considered, effectively, as the 
dose about that point. Qfcoxxrse, the same procedure has. to 


Hans ‘~'2~ ^0 P air 

4 ^ 
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be followed for all points of interest. Finally one can 
get a histogram, say dose vs shield thickness, and from that, 
distribution of doses along the points of interest may easily 
be found. 


Now, the question arise, is, how to choose the 
interval ,pd for a given problem?. As such, there is no specific 
rules or procedures regarding this. Perhaps some experienced 
workers might guess, the procedure to be followed for a problem 
of interest by intution. But in this work, few approaches 
have been made and the results obtained are shown in Figs, 3.1 
to 3,6, 

In the first approach 0.5 and 1.0 were considered for 
^d^ valuGs, i.e. Ad along Z axis (refer Fig. ’2.1 ) in two 
trials. It was found unsatisfactory for even points on the duct 
axis because the values of total dose computed minus the 
unscatter.ed dose computed from equation 3,2 were much higher 
than the unscattered doses (values are not shown here). 

Further some results were obtained for points off 
the axis also using the same approach. Here, the same 2 ^d 
length cannot be considered for every the arc length 



50 


(boundary) is not constant, so it is better to consider equal 
angles, i.e. divide each mfp into eighteen 10 degree; intervals, 
instead of equal heights. Though this helped to remove the 
difficulty faced due to the unequal arc lengths, again, the 
magnitudes of doses were high. 

In the next approach, in addition to^sd„ value, 

Z jC 

i,e, aiong X axis also included to observe the effect on 
dose values. ^d„ = 0.5 and d = 0.5, showed erratic results, 
later when ^ increased from 0,5 to 1.0 (equal to the 
radius of the duct), satisfactory result a have been obtained 
as shown in Table 3,1 for dose points on the duct axis. 

Assuming that the above fixed values of <i5d^ and 
d^ 3 t'e valid for points off the axis also, observation 
have been made for various angles and the results are shown 
in Figs. 3.1, 3.3 and 3.5. For e.g.. Fig. 3.1 shows the 
distribution of dose rates for various shield thicknesses 
at each angular intervals, while Figs, 3.2, 3,4 and 3,6 
show the distribution of dose rates for various points on the 
surface of the shield in each mfps ' . 

The computed multiple scattered doses on the axis 
are shown in Table 3,1 for the three different cases, mentioned 



Table 3,1 Comparis^on of Monte Carlo Dose Rates with 
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in section 3,3^ considered for obtaining these results-. These 
values are compared with the singly scattered doses ccxnputed 
using equat i on ( 3 , 1 ), 

For water at 1.00 Mev photons (survival ratio 1.0 
and for L/Rq = 10; the ratio of maximum shield thickness, 
to duct radius)/ it i s seen that (9) the ratio between 
scattered and unscattered dose is about 40 %^ So, for lead, 
at 1,25 MeV (survival ratio approximately 0.8) and for 1 /Rq 9.0/ 
it is expected that approximately the same 40% will be the 
scattered dose. Table 3,2 shows not much difference between 
the Monte Carlo output and the assumed values. 

As is known that monoeneargetic photons incident 
upon a shield emerges from the shield with a continuous energy 
say from 0 to Eq ( the maximum energy Eq)« So emphasis has 
been shown to count the number of particles falling in an 
energy interval/ say in step of 0.1 Mev from 0,0 to 1,25 Mev, 
while crossing a region or boundary of interest. Further 
this observation v/ill put light on the difference between 
one case and the other (see section 3,3), The observed 
number of particles for different energies at each mfp are 
shown in Figs, 3,7 and 3,8, 
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3.5 DISCUSSIONS 


In the Figs. 3.1, 3.3 and 3.5 curves are not drawn 
for angles 20°, 30° and 40° because the computed values for 
these; angles lie close to each other. The ccmputed value of 
dose at- the fifth mfp at angle 10 is 0.0 in all the three 
cases considered and this point is not shown in all the three 
figures mentioned above. In practice, there is seme contribution 
from a disc soxorce to that point and hence this gives rise 
to doubts about the choice of the Ad^. and D-d^ values chosen. 

A suggested approach is given in Chap. 4. 


Further, number of particles escaping through the 
duct (i.e. part of the uncollided flux) is found to be 
approximately 620 in all the three cases considered for 10^ 
source particles. Since the source energy is only 1.25 MeV 
there is not much difference between curves in Figs. 3.1 and 
3.3, and. Figs. 3.2 and 3.4 because Photo el^txic effect 

which accounts for absorptic^ included. But the 

^ 3 5 and 3.6 which include 

figures: compared with Figs, j.:? 

Photo electric effect show sane difference in the dose values. 

AS mentioned earlier, the number of particles having 




Shield Thickness (mf ps ) 


Rg. 3.7 \fciriation of No. of Particles with Mean 
Free Path. 




No. of Particles 



Fig. 3.8 Variation of No- of Particles with Mean 
Free Path. 
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Table 3.2 Comparision of Monte Carlo Dose Rates 

with modified Dose Rates. (40% above unscattered) 


MFP 

Distance 

In cm 

Dose* 

mR/hr 

Dose M.C. 
case (i) 
mR/hr 

Dose M.C. 
case (ii) 
mR/hr 

Dose M.C. 
case (iii ) 
mJV'hr 

1 

2.547 

mm 

9.392v' 

9.6450 

9.417 

2 

4.0931 

3.9203 

3.770 

3.9740 

3.701 

3 

5.6396 

2.0933 

2.050 

2.1910 

1.986 

4 

7.1861 

1.2967 

1,351 

1.3470 

1.233 

5 

8.7327 

0.8809 

0.892 

0.9033 

0.875' 


* Unscattered Dose + 40% Unscattered Dose 
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particular energy as a function of distance in mfps' are 
shown in Figs. 3,7 and 3,8. Four sets of curves for energies 
0.3, 0.5, 0,6 and 0.9 MeV are shown in those figures. Here 
again, due to absorption in all the four sets^ the inclusion 
of Photo electric effect shows lower: values. The amount 
decreased from that of the other two cases is about 30% to 
65%, At 0.5 MeV,the contribution due to Pair- Product ion is 
more compared to other energies because of the annihilation 


gammas 



CHAPTER- 4 


CONCLUSIONS AND RECOMMENDATIONS 


4.1 CONCLUSIONS 

The following conclusiona can be drawn based upon 
the results of the present study. 

1 • Results have been obtained by using Monte Carlo 

methods, for the investigation of photon transport 
in a chosen shield geometry shown in Fig, 2,1. 

Ccmputed Gamma ray dose rates along the axis of the 
duct are shown in Table 3,1, Those values have good 
agreement with the results obtained by Tsuro et.'ai(ll), 

2. No ccmparision has been made for the doses computed 

at points off the duct axis because neither analytical 
solutions nor experimental results are available, 

3 , In the simulation of photon transport^ all the three 
events Photo electric effect, Pair-Production and 
Compton scattering are included and the effects of 
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absorption due to Photo electric effect and P^ir— 
production are found clearly frcm Figs. 3.7 and 3,8. 

4, From the values obtained for points in the source 

plane, (refer section 3,5), one can conclude that 
the values of a dg and -f-sd^ are not reliable, 

4.2 RECOMMENDATIONS 

1, For the dose calculations, the values of /^d^ andAd^ 
may be changed or the whole scheme of particles 
scoring should be modified suitably to calculate the 
flax accurately. This can be done, for e.g,, as 
follows. Divide the perimeter of the semi circle 
into 18 intervals, try to make a circle about each 
such interval so that the area under each circle 
remains constant and hence, the number of particles 
crossing each circle can be calculated. 

2, This investigation can be carried out by including 
seme biasing techniques for path length selection 
and depending upon the interest, directional biasing 
may also be included. 
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3 , Simulation approach can be extended to shields made 

of hoterogeneous materials or to multilayered shields. 
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